Observational Constraints on Pulsar Wind Theories 



in 
o 
o 



o 



> 
in 
o\ 
in 
in 
o 
in 
o 

Of 

6 

GC 



13 



J.G. Kirk 

Max-Planck-lnstitut fur Kernphysik, Postfach 10 39 80, 69029 Heidelberg, Germany 

j ohn . kirk@mpi -hd . mpg . de 



ABSTRACT 

Two-dimensional, relativistic, MHD simulations of pulsar-wind powered nebulae provide 
strong constraints on the properties of the winds themselves. In particular, they confirm that 
Poynting flux must be converted into particle energy close to or inside the termination shock 
front, emphasising the puzzle known as the cr paradox. To distinguish between the different 
possible resolutions of this paradox, additional observational constraints are required. In this 
paper, I briefly discuss two recents developments in this respect: the modelling of high time- 
resolution optical polarimetry of the Crab pulsar, and the detection of the pulsar/Be star binary 
PSR 1259-63 in TeV energy gamma-rays. 

Subject headings: Stars: pulsars - MHD - Radiation mechanisms 



1. Introduction 

Our view of the nebulae powered by pulsar 
winds has matured dramatically in the past few 
years, largely as a result of high resolution im- 
ages in the X-ray band (Weisskopf et al. 2000; 
Gaensler et al. 2001 ; Helfand et al. 2001 ; Gaensler 
et al. 2002; Lu et al. 2002; Gaensler et al. 2003; 
Pavlov et al. 2003) complemented, in the case of 
the Crab Nebula, by high resolution optical im- 
ages (Hester et al. 1995, 2002). This has moti- 
vated at least three groups of researchers to per- 
form 2D, relativistic MHD simulations (Komis- 
sarov & Lyubarsky 2003, 2004; Del Zanna et al. 
2004; Bogovalov et al. 2005). The pulsar wind 
provides the inner boundary conditions for these 
simulations, which are concerned with the sub- 
sonic flow downstream of the wind's termination 
shock front. Contact with observation is estab- 
lished by computing the synchrotron emissivity of 
the flows, assuming injection of relativistic elec- 
trons with a prescribed spectrum at the termina- 



tion shock. Overall, the observed synchrotron 
morphology, with its rings, asymmetric torii and 
jet structures, is fairly convincingly reproduced, 
although some of the details still evade explana- 
tion (Shibata et al. 2003; Mori et al. 2004), and 
some interesting phenomena demand a kinetic ap- 
proach (Spitkovsky & Arons 2004). There is a 
consensus on the requirements these simulations 
place on the properties of the pulsar wind: the 
energy flux must be concentrated in the pulsar's 
equatorial plane, and, at least in the case of the 
Crab, the fraction of this energy transported as 
Poynting flux must be small. This fraction is mea- 
sured by an average over the surface of the shock 
of the cr parameter, defined as the ratio in the co- 
moving frame of the magnetic enthalpy density 
B 2 /4n to the enthalpy density of the plasma, in- 
cluding its rest-mass. For the Crab Nebula, ob- 
servations suggest cr ~ 10~ 3 . This is puzzling, 
because, on the one hand, standard theories of 
pair production by pulsars (Daugherty & Hard- 
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ing 1982; Hibschman & Arons 2001a,b) imply a 
large value of cr at the point where the wind is 
launched, and, on the other, because a relativistic, 
ideal MHD wind does not appear to permit a sig- 
nificant reduction of cr during propagation (Bo- 
govalov 2001). This puzzle is not new (e.g., Ken- 
nel & Coroniti 1984), and has been termed the "<x 
paradox", or "cr problem" 

Several resolutions of the paradox have been 
proposed. It is possible that magnetic collima- 
tion may indeed result in the desired conversion, 
if the initial conditions in the wind are sufficiently 
anisotropic (Vlahakis 2004). Alternatively, the 
conversion may take place in a thin layer associ- 
ated with the termination shock itself (Lyubarsky 
2003), or may even be to some extent hidden in 
the downstream (nebular) flow (Begelman 1998; 
Shibata et al. 2003). 

A hypothesis that has recently received at- 
tention is that dissipation within the wind zone 
causes the conversion. Coroniti (1990) and 
Michel (1994) originally pointed out that this 
might happen at the boundary between layers of 
magnetic field of opposite polarity in the "striped 
wind". However, the associated acceleration of 
the wind causes a significant dilation of the dis- 
sipation timescale, an effect omitted in the earlier 
calculations. Assuming a gradual acceleration 
process, it is possible to calculate the evolution 
of the entropy wave that constitutes the stripes in 
a short wavelength approximation (Lyubarsky & 
Kirk 2001). Solutions have been found (Kirk & 
Skjseraasen 2003) in which the Lorentz factor T 
varies as a power of the radius r: 

r oc r q with 1/3 < q < 1/2. (1) 

The dissipation proceeds fastest in the solution 
with the slowest acceleration, T oc r 1/3 . Com- 
plete conversion of magnetic into particle energy 
before reaching the termination shock was found 
to be marginally possible in the case of the Crab 
pulsar wind. 

Which, if any, of these possibilities will pro- 
vide the definitive resolution of the paradox can 



only be decided with the help of observational 
constraints. These are difficult to find. The wind 
itself, if mostly cold, does not radiate synchrotron 
radiation and so it is difficult to determine its 
structure. Recently, however, there have been two 
observational developments that do appear to hold 
out the promise of useful constraints. 

In the first, high time-resolution optical po- 
larimetry of the Crab pulsar have become avail- 
able that find a plausible explanation under the hy- 
pothesis that the high energy pulses from the Crab 
pulsar are formed in the pulsar wind. If this inter- 
pretation is correct, the high-energy pulsed radi- 
ation from this and other pulsars may provide us 
with a direct view of the conversion process itself. 

The second development concerns the detec- 
tion in TeV gamma-rays of the binary system 
PSR 1259-63. This offers a diagnostic of the pul- 
sar wind as it enters its termination shock. Be- 
cause the stand-off distance of the shock from 
the pulsar varies with binary phase, it may prove 
possible to extract the dependence of the wind's 
Lorentz factor on distance from the pulsar, thus 
testing the acceleration models. 

2. Polarisation of optical pulses from the 
Crab pulsar 

If the hot particles between the stripes of a rel- 
ativistic, radial wind emit radiation, the effect of 
Doppler boosting will cause it to appear to the ob- 
server to be pulsed, provided that the radius r of 
the radiating zone satisfies 

— < r 2 (2) 

where is the radius at which the corotation 
speed would become luminal and T is the bulk 
Lorentz factor of the wind. This led to the sug- 
gestion (Kirk et al. 2002) that the twin-pulse pro- 
file observed in the Crab pulsar, as well as in the 
gamma-ray emission of other pulsars, is a man- 
ifestation of the fact that two sections of the cur- 
rent sheet are visible per period in the striped wind 
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Fig. 1. — The parameter space of pulsar inclination an- 
gles £ (between the rotation axis and the line of sight) and 
a (between the rotation axis and the magnetic axis). Lines 
joining the lower left with the upper right comers (blue and 
red on-line) are contours of constant phase separation A of 
the twin pulses produced in the striped wind model (equally 
spaced by 0.01 between A = and A = 0.49). Lines join- 
ing the upper left with the lower right corners are lines of 
equal "impact parameter", i.e., minimum angle between the 
line of sight and the magnetic axis. The diagonal (red on- 
line) corresponds to a line of sight that passes directly over 
the magnetic pole, the other two pairs to minimum angles of 
5° (cyan on-line) and 10° (yellow on-line). The spacing of 
the Crab pulses, 0.4, corresponds to the labelled contour (red 
on-line) and the £ at 60° and a » 60° angles favoured by the 
emission model occur at the intersection of this contour with 
the "over pole" diagonal. 



model. If this interpretation is correct, it con- 
strains the the angles f between the rotation axis 
of the pulsar and the line of sight and a, between 
the magnetic and rotation axes, as illustrated in 
Fig. 1. 

In the Crab pulsar, the phase spacing of the the 
high-energy pulses is 0.4 of a period. Combining 
this with the evidence, gathered from the nature of 
the pulsed radio emission (see Kirk et al. 2002), 
that the observer passes almost directly over the 
magnetic pole, fixes both to be roughly 60°. This 



Fig. 2. — Model light curves of the intensity, degree of po- 
larisation and position angle of the pulsed synchrotron emis- 
sion compared to measurements of the Crab pulsar (Kell- 
ner 2002; Kanbach et al. 2003) (from Petri & Kirk (2005)). 
Models with Lorentz factor T = 20 (solid line, red on-line) 
and 50 (dotted line, cyan on-line) are shown. The inclination 
of the line of sight equals the obliquity: a = f = 60° and 
the position angle of the projection of the pulsar's rotation 
axis was set to 124° (Ng & Romani 2004). The bottom pan- 
els show the dependence on phase of the modelled magnetic 
field and particle density in the comoving frame. 



agrees with an independent determination of £ = 
61.3 ±0.1 from the morphology of the X-ray torus 
(Ng & Romani 2004). 

High time-resolution polarimetry of the Crab 
pulsar (Kellner 2002; Kanbach et al. 2003) places 
additional, independent constraints on this model. 
The high-energy pulsed emission (optical to 
gamma-ray) originates as the synchrotron radi- 
ation of relativistic particles accelerated primar- 
ily in the current sheets. The magnetic field re- 
verses direction across a sheet, but, in reality, this 
must occur over a region of finite width. Within 
such a region, the density of hot particles will 
rise, but the magnetic field intensity should de- 
crease. It may become turbulent, and, instead of 
a pure neutral sheet, any residual poloidal field 
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component could cause a sweep of the field direc- 
tion rather than a discontinuous jump. To some 
extent, the population of heated particles must 
also leak into the stripes surrounding the sheet, in 
which case, the toroidal magnetic field component 
should dominate and manifest itself in the polari- 
sation direction. These effects have recently been 
modelled (Petri & Kirk 2005) and compared to the 
data. The results, displayed in Fig. 2 show that it 
is possible to obtain reasonably good agreement 
for the intensity and degree of polarisation. For 
the angle of polarisation, very good agreement is 
achieved with the observed sense and size of the 
polarisation sweeps through the pulses. Particu- 
larly noteworthy, however, is the excellent agree- 
ment between the observed orientation of the lin- 
ear polarisation vector between the pulses and the 
predicted orientation. Within the model, this is 
not sensitive to the choice of parameters, being 
coincident with the projection of the rotation axis 
of the pulsar on the sky. This is determined in- 
dependently of the polarisation observations from 
measurements of the X-ray morphology of the 
nebula (Ng & Romani 2004). 

Thus, the "smoking gun" of magnetic energy 
dissipation in pulsar winds may already have been 
detected (see Arons 2004). However, much work 
will be needed to develop particle acceleration 
models in these sheets before it becomes possi- 
ble to estimate the calibre of the "gun", where it 
fires, and how the trigger mechanism works. 

3. High-energy emission from PSR 1259-63 

The prospect of additional observational con- 
straints on the structure of a pulsar wind has 
been opened up by the recent detection of high- 
energy (TeV) gamma-rays from the Be-star/pulsar 
binary system PSR 1259-63 by the H.E.S.S. array 
of imaging Cerenkov telescopes (Schlenker et al. 
2005; Aharonian, F, et al., HESS Collaboration 
2005). This binary comprises a pulsar of period 
48 ms in an eccentric, 3.4 year orbit about a lu- 
minous Be star. Both stars presumably drive off 
winds, which interact, producing X-rays (Tavani 
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Fig. 3. — Predicted spectrum of inverse Compton scat- 
tered photons from the shocked pulsar wind plasma in the bi- 
nary system PSR 1259-63, 12 days prior to periastron (Kirk 
et al. 1999), and the spectrum observed by the H.E.S.S. array 
of imaging Cerenkov detectors in February 2004 (Schlenker 
et al. 2005; Aharonian, R, et al., HESS Collaboration 2005) 

& Arons 1997). Furthermore, both winds origi- 
nate on rapidly rotating objects and are thought to 
have strong pole to equator variations. Modelling 
of the unpulsed radio emission suggests that the 
axis of rotation of the Be star is not aligned with 
that of the binary (Ball et al. 1999), which compli- 
cates the dependence of the interaction on orbital 
phase. Relativistic electrons from the shocked 
pulsar wind are thought to be responsible for the 
X-rays, in which case they must also upscatter UV 
photons from the Be star to TeV energies. Mod- 
elling of this process led to the prediction shown 
in Fig. 3 for the multi-wavelength spectrum. The 
details of the flow pattern of the emitting plasma 
are unknown, so that it cannot be decided a pri- 
ori, whether the radiating electrons suffer pre- 
dominantly radiative or adiabatic losses. How- 
ever, provided the electron spectrum can be "cal- 
ibrated" using the synchrotron X-rays, this does 
not greatly influence the predicted TeV flux. The 
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relative intensities of the synchrotron and inverse 
Compton components, on the other hand, depend 
sensitively on the magnetic field strength in the 
source. Fig. 3 shows the predicted flux for a mag- 
netic field of 0.3 G, as indicated by models of the 
unpulsed radio emission (Ball et al. 1999), and 
in agreement with estimates from the spin-down 
rate of the pulsar. The remarkable agreement 
in both spectrum and normalisation indicates not 
only that the emission mechanism was probably 
identified correctly, but also that the X-ray inten- 
sity and magnetic field strength vary little from 
one periastron passage to another. 

In addition to the spectrum, the H.E.S.S. ob- 
servations also revealed substantial fluctuations in 
the light curve of this system on a timescale of 
one day. These cannot be accounted for within a 
model such as that of Kirk et al. (1999), which as- 
sumes spherically symmetric stellar winds. How- 
ever, fluctuations are certainly to be expected if ei- 
ther or both winds is asymmetric about the normal 
to the plane of the orbit, especially close to peri- 
astron, where the true anomaly varies rapidly. A 
correlation with the unpulsed radio flux (Johnston 
et al. 2005) would be important in this connection, 
possibly indicating a connection with the passage 
of the pulsar through the equatorial outflow zone 
of the Be star (Kawachi et al. 2004). 

The model underlying the prediction in Fig. 3 
assumes the Lorenz factor of the pulsar wind just 
before it enters the termination shock located be- 
tween the two stars is a few times 10 7 . At peri- 
astron, this point is, at most, a few thousand light 
cylinder radii from the pulsar, and can grow to a 
few hundred thousand at apastron — minute com- 
pared to the 10 9 rL stand-off distance of the termi- 
nation shock of the Crab pulsar. At first sight, 
therefore, the relatively high Lorentz factor ap- 
pears to favour more rapid conversion of magnetic 
energy to particle energy than could be accommo- 
dated in the gradual acceleration models of Kirk 
& Skjseraasen (2003). But the extent to which it 
will be possible to test pulsar wind acceleration 
theories, for example by placing firm limitations 



on the radial variation of the Lorentz factor, will 
depend on the availability of simultaneous X-ray 
and hard gamma-ray observations throughout the 
orbit, as well as the development of more sophis- 
ticated spectral modelling techniques. 
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